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Abstract 

We calculate the sub-gap density of states of a disordered single-channel 
normal metal connected to a superconductor at one end (NS junction) or at 
both ends (SNS junction). The probability distribution of the energy of a 
bound state (Andreev level) is broadened by disorder. In the SNS case the 
two-fold degeneracy of the Andreev levels is removed by disorder leading to 
a splitting in addition to the broadening. The distribution of the splitting is 
given precisely by Wigner's surmise from random-matrix theory. For strong 
disorder the mean density of states is largely unaffected by the proximity to 
the superconductor, because of localization, except in a narrow energy region 
near the Fermi level, where the density of states is suppressed with a log- 
normal tail. 

PACS numbers: 74.80.Fp, 72.15.Rn, 73.63.Rt 
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I. INTRODUCTION 



Several recent works have identified and studied deviations from mean-field theory in 
the sub-gap density of states of a normal metal in contact with a superconductor.EHi The 
excitation spectrum below the gap of the bulk superconductor consists of a coherent su- 
perposition of electron and hole excitations, coupled by Andreev reflectioni at the normal- 
metal-superconductor (NS) interface. The energy of these Andreev levels fluctuates from 
sample to sample, but such mesoscopic fluctuations are ignored in mean-field theory. Be- 
cause of these fluctuations, the ensemble averaged density of states acquires a tail that 
extends below the mean-field gap, vanishing only at the Fermi level (zero excitation energy 
e) . The fluctuations become particularly large if the size of the normal metal is greater than 
the localization length. 

The purpose of this paper is to analyze an extreme case of complete breakdown of 
mean-field theory, which is still sufficiently simple that it can be solved exactly. This is 
the case of single-mode conduction through a disordered normal-metal wire attached to a 
superconductor. The localization length in this geometry is equal to the elastic mean free 
path i, so that the wire crosses over with increasing length L from the ballistic regime directly 
into the localized regime — without an intermediate diffusive regime. Perturbation theory 
is possible in the quasiballistic regime i ^ L, but for i < L an essentially non-perturbative 
approach is required. We will use an approach based on a scaling equation f also known as 
invariant embedding), that has proved its use before in different contexts.l^t 

We will contrast the quasiballistic and localized regimes, as well as the two geometries 
with a single superconducting contact (NS junction) or with two superconducting contacts 
at both ends of the normal metal wire (SNS junction). If we assume that the two supercon- 
ductors have the same phase, so that there is no supercurrent flowing through the normal 
metal, then the Andreev levels of the SNS junction are doubly degenerate in the absence of 
disorder. This degeneracy is broken by disorder. We find that for weak disorder the proba- 
bility distribution of the splitting is given precisely by Wigner's surmise from random-matrix 
theory.0 (The spectra of chaotic systems have spacings described by Gaudin's distribution, 
which is close to, but not identical with Wigner's surmise. 0) 

In the localized regime the fluctuations of the Andreev levels become greater than their 
spacing, and they can no longer be distinguished in the mean density of states, which 
decreases smoothly to zero on approaching the Fermi level. The energy scale for this soft 
gap is exponentially small because of localization, given by Eg = {hvp / t}e~^/^ . The decay 
of {^{s)) for e <^ Eg has a log-normal form oc exp — ^ In^ (e^/e) . Such log-normal tails are 
characteristic of rare fluctuations in the localized regimeEJ and have appeared recently in 
the context of the superconductor proximity effect .i 
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FIG. 1. Geometry of the NS and SNS junctions. 

II. QUASIBALLISTIC REGIME 
A. NS junction 

The NS junction consists of a piece of normal metal of length L connected at one end to 
a superconductor (see Fig. ^a). The width of the normal metal is of the order of the Fermi 
wave length such that there is a single propagating mode at the Fermi energy Ep. We 
assume an ideal junction, without any tunnel barrier and with Ep much greater than the 
superconducting gap Aq. An electron incident on the superconductor with energy e < Aq 
above the Fermi level is then Andreev reflected as a hole at energy e below the Fermi level, 
with the phase shift 

(/)A = - arccos (e/Ao), -7r/2 < 0^ < 0. (2.1) 

We wish to know at which e a bound state (Andreev level) will form in the normal metal. 

The electron and hole components of the wave function ip{x) = {u{x),v{x)) satisfy the 
Bogoliubov-de Gennes (BdG) equation^ 

= eij, (2.2) 

where Tio = —{h^/2m)d'^/dx'^ + V{x) is the Hamiltonian of the normal metal (with disorder 
potential V") and A{x) = Aq6{—x) is the superconducting gap (which vanishes in the normal- 
metal region x > 0). For narrow junctions (width much less than the superconducting 
coherence length = Hvp/Aq) the depletion of A(x) on the superconducting side may be 
neglected, hence the step function 6{—x). At the closed end x = L of the normal metal we 
impose the boundary condition ip{L) = 0. 

In this section we address the quasiballistic regime of mean free path i ^ L. We can 
then treat as a small perturbation on the ballistic bound states 




3 



Hx) = . r;"^^^\v^^^%'"/^^ 1 V 0<x<L, (2.3a) 

y/Z \si'^[{kF - k){x - L) - nn] J ^ ' 

i( \ ^ { sinf/ci.x - (A;f + A;)^] \ { ^ ■ a.\ QK^ 

= • r/ ('z. 7^r 1 ~7~sm0^ , x < 0. (2.3b) 



I sin[/ci?x — (fci? — k')L — vm] y \ 6 







The normalisation constant is Z = L — sin^yi for kpL ^ 1. (We denote kp = mvp/h = 
2n/Xp.) The wave number k = e/Tivp should satisfy the quantization condition 

2fcL + 0^ = 7rn, n = 0,1,2,.... (2.4) 

The total number of Andreev levels within the gap is 2L/7r^o for L ^ ,^o- (There remains 
one level if L <C ^o-) 

To first order in V the energy level is shifted by the matrix element 



5e 



/ dxV{x) u{xf-v{xf . (2.5) 

Jo L -I 



We assume a potential with a short-range correlation, expressed by 

(r(x)) = 0, ly{x)V{x')) = ^5{x-x'), (2.6) 

where (■ ■ ■) stands for the disorder average. It follows that the distribution of an Andreev 
level around its ballistic value is a Gaussian with zero mean, {5e) = 0, and variance 

2£(2L-eo/sin0^)' 

By way of illustration, we show in Fig. ^ the mean density of states of an NS junction 
containing three Andreev levels {^o/L = 0.24) with mean free path i = 12L. The Gaussian 
given by Eq. ( p.7|) agrees very well with the numerical solution of the BdG equation (data 
points). 

We briefly explain the numerical method. The BdG equation is solved numerically on a 
one-dimensional grid (lattice constant a) by replacing the Laplacian by finite differences and 
truncating the Hamiltonian matrix in the superconducting region, where the wave function 
is evanescent for energies in the superconducting gap. The resulting tight-binding model has 
nearest-neighbor coupling 7 = /2'ma? (band width 47). We set Ep = 7 and Aq = 0.1 7, 
corresponding to Ap = 6a and ^0 = IOa/So. The disorder is modelled by a random on-site 
potential which is uniformly distributed in the interval (— VF, W). The mean free path from 
the Born approximation, I = SEplAj — Ep)a/W'^, was found to fit well to the prediction 
of one-dimensional scaling theory for the mean inverse transmission probability, (T~^) = 
|[1 -|-exp(2L//)], in the complete range from the quasiballistic to the localized regime. (The 
localization length ^ is related to the mean free path by ^ = 2£, cf. Ref. ^) This allows for a 
parameter-free comparison of the analytical and numerical results for the ensemble-averaged 
density of states. 
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FIG. 2. Mean density of states (in units of i^q = IL/irhvp) of a quasiballistic NS junction. The 
Gaussian with variance given by Eq. ( |2.7D (sohd curves) is compared to the numerical solution of 
the BdG equation (data points). 



B. SNS junction 

The quasiballistic regime in an SNS junction (Fig. |l|b) is qualitatively different from the 
NS case of the previous subsection. The reason is the double degeneracy of the unperturbed 
Andreev levels. This degeneracy exists if the phase of the order parameter in the two 
superconductors is the same, which is what we assume in this paper. Let us examine the 
splitting of the Andreev levels by the disorder potential. 

The SNS junction has energy gap 

A{x) = AoOi-x) + AoOix - L). (2.8) 

The quantization condition reads 

kL + (f)A = n7i, n = 0,1,2,.... (2.9) 

There are L/n^o Andreev levels (for L ^ ^q), each level being doubly degenerate. We choose 
the two independent eigenf unctions ipii^) such that they carry zero current. They are given 

by 

= ( cos'rix-"^ L) ) ("i ^ °' ^^'^^""^ 

Mx) = 4= ( "uf^'t^ 1 V 0<x<L, (2.10b) 

^/Z' \cos[{kF - k)x - 4)a] J ^ ' 

, . I ( cosikpx + kL)\ (x-L \ (o ^(^ \ 

%pAx) = —=\ }, , ; exp — - — sm0A , x > L, (2. 10c) 
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FIG. 3. Distribution of the splitting s of the first pair of Andreev levels in an SNS junction 
with ^q/L = 0.24. The solid curves are our theoretical expectation from Eq. ( p. 11 ), the data points 



result from the numerical solution of the BdG equation. The inset shows the numerical data for 
the mean density of states. 

and ip-lx) is obtained by replacing cosine by sine. The normalization constant is now 
Z' = L- ^o/ sin0^. 

To first order in V the levels are splitted symmetrically around the ballistic value, by 
an amount ±|s. The basis ( p.lO|) is chosen in such a way that the off-diagonal elements of 



the perturbation vanish. The shift of each level can then be calculated from Eq. ( |2.5| ) using 
the corresponding eigenf unction. We again calculate the probability distribution P{s) of the 
level splitting using Eq. (|2.6| ). The result is 



US I TTS^ 



with average splitting 



(s) = A,/^^5^^1^±M^. (2.12) 



L - ^o/ sin 0A 

We recognize Eq. ( |2.11|) as Wigner's surmise of random matrix theory.0 



In Fig. ^ we compare Eq. ( p.ll| ) with numerical data. The agreement is excellent for 
a range of mean free paths in the quasiballistic regime. The mean position of the splitted 
levels fluctuates only to higher orders in Lji. This makes it possible to resolve the splitting 
in the mean density of states (see inset in Fig. 
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III. LOCALIZED REGIME 



A. NS junction 

In order to go beyond the quasiballistic regime into the locahzed regime L ^ i we write 
the quantization condition for the Andreev levels in an NS junction in the form 

r{e)r{-e)*e^''^^ = 1, (3.1) 

where the reflection amplitude of the disordered normal metal. [The hole has 

reflection amplitude r*{—e).] In terms of the phase shifts we have 

^ m ~4>{-e) ^ ^^^^^ = vrn, n = 0, 1, . . . . (3.2) 
The density of states i'{e) = Y.n^{^ ~ ^n) is related to the scattering phase shifts bylii 

z/(£) = --^Im lnsin$(e + iO+), (3.3) 
vr ae 

where O"*" denotes a positive infinitesimal. The imaginary part of the logarithm jumps by vr 
whenever sin $(e) changes sign, hence it counts the number of levels below e. The derivative 
with respect to e then gives the density of states. It is convenient to write Eq. ( |3.3| ) as a 
Taylor series, 



, , I d . ^ 1 



, , $ + Im y Ae^^'"* , (3.4) 
7rde\ ^1 m J 



which converges because $(£ + iO"*") is equivalent to ^{s) + iO~^. 

We seek the disorder averaged density of states (i^(e)). One way to proceed is by means 



iique 



of the Berezinskii technique.lijtJ An alternative way, that we will follow here, is to start 



from the scaling equatioroQ for the probability distribution P{(Pn) of the phase shift 0jv = 
i [(/)(£:) — (f>{—e)]. This equation has the form 

dP d f 2e 1 d ^ \ ^ , , 

The initial condition is lim^^o -P(0Af) = 5(0Af)- 

The first moment satisfies d{(j)N)/dL = 2e/hvF, hence 

(^.) = — . (3.6) 

Multiplication of Eq. ( p.5| ) by expi2im0Ar) and integration over 0Ar from to tt yields a set 
of recursive differential equationsllj for the moments Rra = (e^*™''^^), 

= — (-Rm+l + — 2i?m) + -^^mRm, (3.7) 
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FIG. 4. Mean density of states of an NS junction from the quasiballistic into the locahzed 
regime. The sohd curves have been computed from Eqs. ( p.4| ) and (3.7). The dashed curves are a 
numerical simulation of random disorder in the BdG equation. 



with the initial condition -Rm(O) = 1. We solve this set of equations by truncating the vector 
i?25 • • • Rm) at a sufficiently large value of M ^ 400 and diagonalizing the corresponding 
tri-diagonal matrix. From Eq. ( |3.4]) we then find the mean density of states. 

The result is shown in Fig. | for C,o/L = 0.24 and ratios £/L ranging from the quasiballistic 
regime to the localized regime. Agreement with the numerical solution of the BdG equation 
is excellent over the whole range. 

In the localized regime L ^ i the individual Andreev levels can no longer be distinguished 
in the mean density of states, because the broadening of the levels becomes greater than 
the spacing. In this regime we distinguish two energy ranges, e ^ Eg and e Eg, where 

For energies higher than Eg we may use the L —>■ oo limit of the distribution P{(j)N), 
obtained by setting the left-hand-side of Eq. equal to zero. The resulting moments are 



lim Rr, 



da e 



a 



a 



iEi 



(3.8) 



We then calculate the mean density of states from Eq. 



2L 



+ 



1 



, with the result 



a — iuj o"(l — 6^**^-4) — zcj 



(3.9) 
(3.10) 



The first term on the right-hand-side of Eq. (|3.9|) is the energy independent density of states 
uq in an isolated normal metal. The main effect of the superconductor for e ^ Eg is an 
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enhancement of the density of states close to the gap Aq of the bulk superconductor (second 
term). The third term is negative for sufficiently small e and is a precursor of the soft gap 
near the Fermi level. For ^ and e <^ Tivpjl the reduction term j[e) can be simplified 
as 



In- 



lei 
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(3.11) 



where 7 ~ 0.58 is Euler's constant. 

Near the Fermi level, for e <^ Sg, the mean density of states vanishes as a result of the 
proximity to the superconductor. This "soft gap" appears no matter how strongly localized 
the normal metal is. The coefficients Rm may now be treated as analytical functions of the 
parameter 

R^ = R{z). (3.12) 



Taking the limit e 



we deduce from Eq. 



.d^R 



3D the partial differential equation 



zR, 



(3.13) 



^dR ^ 

dL dz^ 

with initial condition lim^^^o -^(-2) = 1- This difi'erential equation has been studied before 
in the theory of one-dimensional localizationjlllll but not in connection with the proximity 
effect. The result for the mean density of states, derived in the Appendix, is given by 

2i 



ui / ui 

m — 

2L \ 2L 



(3.14) 



Innhvp/ei = InnEg/e + L/i. The leading logarithmic asymptotic of this expres- 



where u 

sion in the limit e Eg has the log-normal tail 



{i^{s)) oc exp 



4L e 



9- 



(3.15) 



The same log-normal tail was found in Ref. § for a many-channel diffusive conductor. In 
that case the factor i/L is replaced by the Drude conductance of the normal metal and the 
energy scale Eg is replaced by the Thouless energy hD/L^ (with D the diffusion constant). 
In our single-channel localized conductor nether the Drude conductance nor the Thouless 
energy play a role. 



B. SNS junction 

In contrast to the quasiballistic regime, the NS and SNS junctions are similar in the 
localized regime. (At least for the case of zero current through the SNS junction considered 
here.) Unfortunately, there exists no simple scaling equation as Eq. (|3.5| ) that can describe 
the density of states of the SNS junction. We therefore rely on the numerical solution of 
the BdG equation. In Fig. |^ we show that the mean density of states of an NS junction of 
length L is close to that of an SNS junction of length 2L. This factor of 2 has an obvious 
explanation in the ballistic regime [compare Eqs. ( |2.4| ) and ( p.9| )], but it is remarkable that 
it still applies to the localized regime. 
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FIG. 5. Numerical calculation of the mean density of states of an NS junction (solid) and SNS 
junction (dashed) in the nearly localized regime. The length of the SNS junction is twice that of 
the NS junction. (The weak oscillations are remnants of Andreev levels, that will disappear \i Ljl 
is increased further.) 



In summary, we have calculated the effect of disorder on the spectrum of Andreev levels 
in single-channel NS and SNS junctions. The non-perturbative effects of localization in the 
one-dimensional case can be studied exactly, at least in the NS geometry. Our research is of 
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APPENDIX A: DERIVATION OF THE LOG-NORMAL TAIL 



The differential operator on the right-hand-side of Eq. ( ^.131 ) has eigenfunctions 



(Al) 
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where Kp{z) is the modified Bessel function, such that 



,9^ 



p 



(A2) 



The solution to Eq. ( |3.13| ) with the initial condition \im.L^Q R{z) = 1 is 

V sinh('7ri//2) 



R{z) = h{z) + j du 



7r(z/2 + 1) 



(A3) 



To obtain the density of states of the NS junction it is convenient to define the inverse 
Laplace transform 



200+0+ 



-joo+0+ 



(A4) 



From Eq. ( p.4|) we find for e <^ Eg the mean density of states in terms of the function F, 



He)) 



' si 

TchVi 



(A5) 



Our aim is to find the asymptotic form of F{X) in the limit A — > 0. The inverse Laplace 
transform of the modified Bessel functions in Eq. ([A 3D can be found in Ref. IT^. We obtain 



oo 



X ■ 



iFi(| + f,l + 2z/, -4A 
2VvF(l -iz/)r(ii//2) 



(A6) 



where Fq{X) = exp(— 4A). The integrand has a single pole u = —i in the lower half of the 
complex plane and the residue from this pole cancels the term Fq. Let us shift the contour 
by the transformation (i£/L) In (1/A) and consider the limit A <^ e~^^^. In this 

limit the contour is shifted through the pole so that the term Fq is cancelled. Moreover, the 
hypergeometric function iFi can be replaced by unit in this limit. Thus, we end up with 
the integral 



F(A) 



X 



2v^ 



2 oo 



L 



InA r 



due 
iu 



2L 



InA 



-1 



(A7) 



The asymptotic form of this integral in the limit A -C e~^^^ can be found by evaluation of 
the expression in square brackets in the point = and calculation of the Gaussian integral. 
Using the asymptotic formula for the Euler gamma function one obtains the mean density 
of states given in Eq. ( p. 141 ). 
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